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Abstract 

We describe CIJETl.O, a Fortran program that aiming for the calculation of single- inclusive jet or 
dijet production cross sections induced by quark contact interactions from new physics at hadron 
colliders, up to next-to-leading order in QCD. It covers various contact interactions with different 
chiral and color structures. The code is designed in a way that could be used for fast calculations 
with arbitrary parton distribution functions based on interpolations of the QCD coupling constant 
and parton distribution functions. 
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1. Introduction 

Jet production at hadron colliders provides an excellent opportunity to test perturbative QCD 
(PQCD) and to search for possible new physics (NP) beyond the Standard Model (SM) over a wide 
range of energy scales. Invariant mass distributions of the dijets dijet angular distributions 
inclusive jet px distribution [HI, and other jet observables at the LHC ]6[ have already extended 
current searches for quark compositeness, excited quarks, and other new particle resonances toward 
the highest energies attainable. Among all these measurements, the inclusive jet pr distribution 
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and dijet angular distribution show a great sensitivity to possible quark contact interactions (CIs) 
induced by new physics models. In the SM, Quantum Chromodynamics (QCD) predicts the jets are 
preferably produced in large rapidity region, via small angle scattering in i-channel processes. Also 
the production rates fall off rapidly in large px or large dijet invariant mass region. On the contrary, 
the jet production induced by quark contact interactions is expected to be much more isotropic and 
fall off slower as the increasing of jet px or dijet invariant mass, and thus the distributions at the LHC 
could be largely modified. 

The measurement of quark contact interactions has been used to set limits on the quark composite 
models which have been studied extensively in the literature 0, Q • It is assumed that quarks are 
composed of more fundamental particles with new strong interactions at a compositeness scale A, 
much greater than the quark mass scales. At the energy well below A, quark contact interactions 
are induced by the underlying strong dynamics, and yield observable signals at hadron colliders. 
The newest bounds of A at the 95% confidence level (C.L.) from the CMS collaboration are around 
lOTeV based on 2.2 fb~^ collected data of the dijet angular measurement and 5fb^^ data of the 
single-inclusive jet measurement [^. Previous limits from the Tevatron and LHC can be found in 
Refs. BiiQ. 



In our earlier studies we presented the theoretical calculations of the next-to- leading (NLO) 
order QCD corrections to the dijet production induced by quark CIs, which are used by the exper- 
imentalists to set constraints on the quark compositeness scales In this work we summarize the 
numerical program used for the dijet calculation with extensions for the calculation of single-inclusive 
jet production, which will be used in the upcoming measurement One common issue of the ex- 
perimental analyses at the hadron colliders is that, they require massive repetitive calculations using 
different parton distribution functions (PDFs) to estimate the theoretical uncertainties. And fast cal- 
culation algorithms for arbitrary PDFs, like the implementations of FastNLO [l^l and APPLGRID [l3| 
are highly desirable. Thus in our program we further provide a fast interpolation modular that allows 
for the calculations of using arbitrary PDFs within seconds while maintaining the desirable numerical 
accuracy. 

This document is structured as follows: Section [5] briefly summarizes the theoretical model and 
calculations. Section [3] describes the fast interpolation algorithm used. And Section |4]|6] show the 
inputs, running, and outputs of the program. Section [7] shows some sample results from the program 
and Section m is a conclusion. 



2. Jet cross sections and quark contact interactions 

We consider a subset of quark contact interactions that are the products of clcctroweak isoscalar 
quark currents which are assumed to be flavor-symmetric to avoid large flavor-changing neutral-current 
interactions Q . The effective Lagrangian can be written as 

1 ^ 

^NP^ ^^c,0^, (1) 

1=1 

where A is the new physics scale, Ci are Wilson coefhcients. And the operators Oi in chiral basis are 
given by 

/ 3 3 \ 

01 = SijSkl E (iLcilfiqLcj E qLdkl^QLdl , 

\c=l d=l I 

02 = T^tjT^ti E qLc^l^,qLcJ E Udkl^qbdi 1 , 

\c=l d=l ) 

/ 3 3 \ 

03 qLciit^qLcj E qndki^qRdi , 

\c=l d=l ) 

04 = TfjTli E qLctif,qLcj E qndki^qRdi , 

\c=l d=l ) 
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Os = S.ijSki ^ qRciliiqRcj ^ qRdkl^qRdi 



\c=l d=l / 

Oe = T^jTJJ, ^ qRcilf^qRcj qRdkYqRdi j , (2) 

\C=1 £i=l / 

in which c, d are generation indices and i. j. k, I, a are color indices, and T" are the GcU-Mann 
matrices with the normaUzation Tr(T''T'') = 5°-^/2. Beside of the quark compositeness, the above 
interactions can also arise from various kinds of new physics models, induced by the exchange of new 
heavy resonances, such as Z' models and extra dimensions models [l5l |. 

The jet cross sections (up to NLO) in each kinematic bin, depending on the scale A and Wilson 
coefhcients Ci at input scale A, can be expressed as 



1=1 



1=1 

(AiAi+i(&ii+i + aji+ir))/A^ 

= 1,3,5 

^ {X^Xi{b,i + aur))/A\ (3) 



1=1,2,5,6 



with Xi defined by Ci ~ AnXi, r ~ hi{A/ ^q), and fiQ is an arbitrary reference scale chosen according to 
the kinematic range of the bin (not the same as QCD scales). All a's vanish at leading order (LO). 
Since QCD interaction preserves parity, we should have 6(a)i_2 = fe(a)5,6, ^(1)11,22,12 = ^(0)55, 66, 56, 
and fe(a)i4,24 = 6(a)54_64- More details of the theoretical calculations can be found in fllj . 



3. Grid interpolation 

In general, dependence of the jet cross sections on the PDF and the QCD coupling constant as 
can be expressed as 

i,3 P 

where i,j are the parton flavors, p indicates the order of as, and ^if[r) ^ire the factorization and 
renormalization scales. We choose a 3-dimensional grid of — 2:2 — /i/ with x x Nq grid 
points, which are uniformly spaced in x\, x\ and Inln/xy/Qo with 5 = 0.3 and Qo = 0-3 GeV. We 
choose Nx = 30 and Nq = 15 which are found to be good enough with an interpolation accuracy of 
about 0.1%. Both ends of each dimension are determined according to the kinematic range of the bin 
considered. We use a fourth-order interpolation and rewrite the cross section as 

(Jbm^Y ^ 4,p(ni,?T-2,n3,0(a^($/U/)$fe)|„i,„2,„3, (5) 

p,k ni,ri2,n3 

where ^ ~ fir/fJ-f, and Ui are indices of grid points. The sum is performed over the perturbative 
expansion order of QCD coupling constant (p) and the parton luminosities (k). The parton luminosities 
involved here for up to NLO calculations are 

ij=qq' M' 

$2 = Y ^1^2ft/hi{xi,flf)fj/h2{x2,fJ.f), 

*3 = X! ^l^y^/hiixi, fJ-f)fj/h2{x2, fJ-f) + {hi ^ h2), 
ij=qq' 

$4= ^ xlxlf,/f,,{xi,flf)fj/h2{x2,fJ.f) + {hi^h2), 

i3=qq 
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ij=gq,m 



(6) 



Note that the gg luminosity does not contribute to the cross sections induced by quark contact 
interactions up to NLO. In the code the grid coefficients ^^^2 5 "-Si C) with dependence on ^ 

solved analytically, are calculated and stored in a table file for each kinematic bin, which allows fast 
interpolations of Ubin for arbitrary PDF and rcnormalization scale choices. The grid interpolation we 
used here is based on the original approaches introduced by FastNLO 12 1 and APPLGRID [I^, but 
with slightly modifications to adapt for our calculation. 



4. Program inputs 

The input parameters are specified in the files with the extension . card in the subdirectory data. 
Each line contains a record for an input variable: a character tag with the name of the variable, 
followed by the variable's value and descriptions. 

proinput . card contains various switches that 

• pdf is the name of the PDF file read from LHAPDF used in the calculation, e.g., CTIO . LHgrid 
for CTIO NLO PDFs pLZl]. Since the code has a fast PDF interpolation interface in the final 
step, user can also change the PDF choice there. 

• pdfmember specifies the PDF member used in the calculation as appearing in LHAPDF, e.g., 
for central PDF. 

• ang specifies the two observables of the kinematic bin. Basically the program calculates the 
total cross section in two-dimentional kinematic bins. For aiig=l-3, the first observable of the 
bin is always the dijet invariant mass. The second observable is x = oxpdyi — y2\) for ang=l, 
V* = \Vi ~ 2/2I/2 for ang=2, and y^ax = niax(|yi|, \y2\) for ang=3. While ang=4 is for single- 
inclusive jet production with the first observable as pT of the individual jet and second as 
absolute rapidity Ij/I of the individual jet. 

• mode=f itll/f itcs/f ital controls the calculation scheme. The code first calculates the coeffi- 
cients (a, h) defined in Eq. ^ at both LO and NLO. The differences between the three choices 
are that f itll assumes only left-handed color-singlet operator at input scale (coupling ci). And 
f itcs includes all the color-singlet operators at input scale (couplings ci, C3, C5). f ital con- 
siders all the six operators in Eq. ([T|). User can choose the calculation scheme according to the 
model studied since the most general case mode=f ital consumes much more CPU times. 

• pseed specifies the random-number seed used by the Monte Carlo integration. 

• many is the number of Monte Carlo sampling points specified for the calculation. User may need 
to adjust it according to the numerical accuracy required. For an accuracy of a few per mille, 
recommendation is 3,000,000 or more. 

• ppcollider represents the type of the collider. for pp machine and 1 for pp machine. 

• Sqrts gives the centcr-of-mass energy y/s of the collider in GeV. 

• scalescheme specifies the choice of the hard momentum scale that defines the central value of 
the factorization and renormalization scales. For dijet production (aiig=l-3), corresponds to 
choose the central value as average pr of the two leading jets, and 1 corresponds to choose the 
central value as exp(0.32/*) times the harder jet pr in the two leading jets. While for single- 
inclusive jet production, indicates using the individual jet pT as the central value, and 1 for 
using maximum jet pT in the rapidity region as the central value. 

kininput . dat contains parameters for clustering and selection of jets, in the same format as in 
proinput . card. 
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• j etscheme specifies the je t alg orithm used in the calculations: 1 for the anti-fcy jet algorithm [l 



2 for the cone algorithms [19|, with Rsep=l .3 for modified Snowmass, and Rsep=2 for midpoint 



algorithm. Note here for the NLO calculation the anti-fcr algorithm is equivalent to the kx [20| 
or Cambridge- Aachen [2lj algorithm. 

• recscheme sets the recombination scheme (22j used in jet clustering: 2 is for the 4-momentum 
scheme (recommended for LHC), and 1 for the Et scheme. 

• Rcone sets the distance parameter or cone size for the jet algorithm used. 

• Rsep is only active for cone algorithm as explained above. 

• ycut and ptcut specify conditions for jet acceptance, i.e., the maximal rapidity and lowest pT 
(in GcV) for a cluster to be considered as a jet. 

• yb, ys specify the upper limit on ytoost — \yi + 2/2 1/2 and y* of the dijet system, respectively, 
and are not used for calculations of single-inclusive jet production. 

bininput . card sets up a scan of the two-dimensional bin. Under tagged line massbin, user can 
specify one invariant mass (or pt) bin per line, containing the lower and upper limits (in GeV). And 
under tagged line rapbin, user can write one angular bin per line (corresponding to Xi V* i Vmax, or 
\y\ depending on ang). The code will do a scan over all the two-dimensional bins specified. 

5. Compilation and running 

After unzip the source file, go to the main directory and run make. The executables will be 
generated and moved into directory data. The program uses the VEGAS subroutine from CUBA 
library [2^ for the Monte Carlo integration. 

There are two main executables. dijets4ci_sin is for the calculation using a single CPU core. 
And dijets4cijnul is for the calculation using multicore which is much faster. Before runing 
dijets4cijnul, user needs to first set the number of cores used by running source setcores.sh 
$1, where $1 is the number. Otherwise dijets4ci_mul will use all the cores available. To run the 
code simply by entering . /dijets4ci_sin $tag or . /dijets4cijnul $tag under the data directory, 
where $tag specifies the name of the job. Note that during the run the code needs to write into some 
auxiliary files for the parellel calculation. Also users are suggested to use reasonable amount of cores 
since it costs additional times for the master process to combine results from each cores. Besides, there 
are another two executables in fastCI directory, i.e., ciconv and cixsec, which will be explained 
later. 

6. Program outputs 

At the end of each run, the program generates an output file $tag_f itresults . dat containing the 
calculated coefficients in Eq. ([3]) for the specified PDF, 6 groups in total for each bin (LO and NLO 
with 3 scale choices). Note that the value of /xq for each bin are also written into the output file. 

At the same time, the gird files (one for each bin) are generated and stored in the f astCI/f grid 
directory. Under fastCI directory, user can first run ./ciconv $1 $2 $3 $4 to calculate the co- 
efficients for arbitrary PDFs and scale choices, where argument $1 is the PDF file name, $2 is the 
number of PDF member, $3 is the grid file name as in f grid directory, and $4 is the specified file 
name for the output. The calculated coefficients, 18 groups for each bin (LO and NLO with 9 scale 
choices), are written into file $4. Using the coefficients, as well as /ig value in $4, user can calculate 
the cross sections for arbitrary scale A and Wilson coefficients q according to Eq. Or user can 
run ./cixsec $1 $2 to get the cross sections for A and Aj values specified in file input, where $1 is 
the output from ciconv, and the cross sections are written into $2. 
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7. Sample results 



In this section wc present some representative sample results from the program for both single- 
inclusive jet and dijet production at LHC 8 TeV. For simplicity we just show results for a single model 
with A = 7 TeV and vector-like color-singlet coupling, Ai_5 — —1, A3 = —2, A2,4.6 = 0. We use the 
anti-fc^ jet algorithm with i? = 0.7 for single-inclusive jet production and K = 0.6 for dijet production. 
The scalescheme is set to for single-inclusive jet production and 1 for dijet production. And we 
only present the differential cross sections in the central rapidity region, i.e., \y\ < 0.5 or \y* \ < 0.5 for 
single-inclusive jet and dijet production, which is preferred by the contributions from CIs. In Fig. [1] 
we show comparisons of the LO and NLO jet cross sections. The NLO QCD corrections are found 
to reduce the cross sections in most of the cases except for the extremely large px or large invariant 
mass region. Also in Fig. [1] we plot the ratios of the cross sections that directly from the numerical 
integration and from latter interpolations. We can see very good agreements between them, and 
the interpolation only introduces numerical errors of the order 0.1%, which are negligible in current 
analyses. 

Using the fast interpolation we can obtain the results for arbitrary scale and PDF choices without 
redoing the numerical integration, which is the most time consuming part. In Fig. [2] we plot the scale 
variation envelope of the jet cross sections by varying the factorization and renormalization scales 
independently around the central value by a factor of two. The cross sections are normalized to the LO 
or NLO predictions using the central scale value. We can see that the NLO corrections greatly reduce 
the scale uncertainties, especially the factorization scale uncertainties, for both cases, which makes 
the theoretical predictions more reliable. Furthermore, in Fig. [3] we compare the NLO predictions 
with 68% C.L. PDF errors from CTIO [13, MSTW2008 0, and NNPDF2.3 (H] NLO PDFs. CTIO 
gives higher predictions in both cases, and MSTW2008 shows much smaller PDF uncertainties in the 
extremely large px and large invariant mass region. 
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Figure 1: Jet cross sections induced by the CIs at both LO and NLO, using the central scale choice and CTIO NLO 
PDF. Ratios between the results directly from the numerical integration and the results from latter interpolation using 
the weight grids are also shown. Left and right panel correspond to single-inclusive jet and dijet production respectively. 



8. Conclusions 

In conclusion, this document describes the CIJETl.O program that aiming for the calculation 
of single-inclusive jet or dijet production cross sections induced by quark contact interactions from 
new physics at hadron colliders, up to NLO in QCD. It covers various contact interactions with 
different chiral and color structures. The program allows for parellel calculations of the numerical 
integrations using the CUBA library. Moreover, the program includes a fast interpolation modular 
that could be used for calculations with arbitrary parton distribution functions without redoing the 
time consuming numerical integrations. The program could be used for the experimental analyses of 
the quark compositeness at the LHC. 
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Figure 2: Scale uncertainty envelope of the jet cross sections induced by the CIs at both LO and NLO normalized to 
the central scale predictions, using the CTIO NLO PDF. Left and right panel correspond to single-inclusive jet and dijet 
production respectively. 




Figure 3; PDF uncertainties of the jet cross sections induced by the CIs at NLO normalized to the CTIO central 
predictions, using the central scale choice. Results are shown for CTIO, MSTW2008 and NNPDF2.3 NLO PDFs. Left 
and right panel correspond to single-inclusive jet and dijet production respectively. 
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